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1. Introduction 

Understanding the phase structure of Quantum Chromodynamics (QCD) at nonzero tempera- 
ture (T) and quark chemical potential (/x) is an important subject relevant to many areas of physics, 
including relativistic heavy ion collisions, the early universe, and the interior of neutron stars Ep. 
Although QCD at T > and pL = has been extensively investigated by both analytical methods 
and first-principle lattice simulations, QCD at /I / is much less understood, as the fermion sign 
problem makes Monte Carlo simulations extremely difficult. One exception is QCD in the high- 
density limit (jj. 3> Aqcd) where the QCD coupling is weak: A color superconductor called color- 
flavor locked (CFL) phase was shown to be realized ^ in which color and flavor symmetries lock 
each other and chiral symmetry is spontaneously broken. 

At jU = an intimate relation has been established between the spectral properties of the 
Dirac operator and the spontaneous breaking of chiral symmetry. In 1992 Leutwyler and Smilga 
[Q] succeeded in deriving details of the Dirac eigenvalue distribution from the low-energy effective 
theory of QCD in a finite volume (the £ -regime) and found that the distribution of Dirac eigenvalues 
near zero is governed by the nonvanishing chiral condensate. Soon after their work, Verbaarschot 
et al. discovered that chiral random matrix theory (ChRMT), which possesses the same global 
symmetries as QCD, provides rich information on spectral correlations of the Dirac operator ^ |6|], 
leading to the realization that the way in which the thermodynamic limit of the spectral density near 
zero is approached is universal. This feature enables us to determine the magnitude of the chiral 
condensate in the QCD vacuum with considerable accuracy by matching the numerical results from 
lattice QCD simulations against the exact analytical results obtained in ChRMT [f7|]. In the past few 
years considerable progress has also been made in ChRMT at pL / (see [||] for a review), but the 
high-density region /I 3> Aqcd is largely unexplored so far. 

Recently two of us pointed out a new £-regime specific for the BCS-state of QCD at high 
density, where exact analytical results, including Leutwyler-Smilga-type spectral sum rules char- 
acterized by the fermion gap A, can be derived This explicitly shows that the Dirac spectrum 
at high density is governed by the BCS gap but not by the chiral condensate. In this report we 
apply these lines of analysis to two-color QCD with an even number of flavors. Our principal 
motivation comes from the fact that two-color QCD can be simulated on the lattice even at /I ^ 0. 
Although this property makes the theory a very attractive testing ground for methods and concepts 
developed in the studies of three-color dense quark matter, the related works have focused on the 
Bose-Einstein condensate of the diquarks (see, e.g., [|I0|]). On the other hand, in this report we 
will concentrate on the BCS superfluid, which is the genuine two-color counterpart of the color 
superconductivity (e.g., CFL phase) in the three-color case. Testing our results on the lattice will 
provide the first signature of BCS pairing as well as determine the value of the gap. 

2. Low-energy effective theory 

Let us first construct the low-energy effective theory for dense two-color QCD [fTlj]. We first 
define our notation. The fermionic part of the Lagrangian in Euclidean space reads + 
with the /I -dependent Dirac operator 

= 7 v D v + 7oAi (2.1) 
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and the mass term 

^ = \(l + j5)M + ^{\-y 5 )M^ . (2.2) 

Here, y is a short-hand notation for Nf flavors of two-color Dirac spinor fields transforming in 
the fundamental representation of SU(2) co ] or . The y v are hermitian /-matrices. The covariant 
derivative D v is an anti-hermitian operator so that the eigenvalues of y v D v are purely imaginary. 
M is the Nf x Nf quark mass matrix. For M = and jj, = 0, the fermionic part of the Lagrangian is 
symmetric under \J(2Nf). For pL ^ 0, this symmetry is broken explicitly to SJJ(N/)l X SU(Nf)g x 
U(1)bxU(1)a- A remarkable property of two-color QCD is that the fermion sign problem is 
absent at nonzero /x: Because of the pseudo-reality of SU(2) we have £F(/i)T 2 Cy 5 = T 2 Cy 5 £F(/l)*, 
where C is the charge conjugation operator and T 2 is a generator of SU(2) co i or . Together with chiral 
symmetry, {75, = 0, it follows that if X is one of the eigenvalues of so are —A, A*, 

and —A*. Consequently, the fermion determinant is real and non-negative in two-color QCD with 
an even number of flavors with pairwise degenerate quark masses. 

At n S> Aqcd, perturbative one-gluon exchange indicates that the color antisymmetric channel 
1 (coming from 2 (g) 2 = 3 © 1) is attractive, which implies that the Fermi surface becomes unstable 
and subject to the formation of Cooper pairs, as known in the BCS theory. As a result a gap 
A appears in the spectrum of quasiquarks near the Fermi surface. The Pauli principle forces the 
condensation to occur in the flavor-antisymmetric channel: 

+ ( w) = (s ah (YlyCY 5 I ij Yl) , (2-3) 

where a, b G {1,2} and i,j G {1, . . . ,Nf} are color and flavor indices, respectively. The Nf x Nf 
symplectic matrix / is defined as 




(2.4) 

where 1 is the (Nf/2) x (Nf/2) unit matrix. Hereafter we always assume that Nf is even. 



Some comments are in order. Although the diquark pairing (2.3) has the same quantum num 



bers as that in [10], their physical meanings are different: The pairing (2.3) is a weakly-coupled 



BCS-type condensate whereas that in [ |10| ] is a strongly-coupled Bose-Einstein condensate (BEC), 
and there could be a smooth crossover between the two. The crossover from low to high densities, 
if realized, is a typical BEC-BCS crossover known in condensed matter physics. It parallels the 



idea of quark-hadron continuity [ |12| ] in real (three-color) QCD, which may be explicitly realized by 
the effect of the axial anomaly [O]. However, for two-color QCD with even Nf, the axial anomaly 
never acts as an external field for the chiral condensate, and a new mechanism is necessary to 
account for the crossover phenomenon. 



The condensation ( |2.3[ ) breaks chiral symmetry spontaneously as follows: 

SU(^V / ) L x SU(^V / ) R x U(1) B x U(1) A -» Sp(^V / ) L x Sp(N f ) R . (2.5) 

We used the fact that the U(l ^-anomaly is suppressed at high density. It is a characteristic of 
two-color QCD that Cooper pairs are color singlets and preserve gauge symmetry; the system is 
in the superfluid phase but not in the superconducting phase. We sttess that ( |2.5[ ) is different from 
the way chiral symmetry is spontaneously broken at /I = Jl0[]: SU(27V/) — > Sp(2Nf). (Thus the 
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resulting effective chiral Lagrangian is entirely different as well.) The Nambu-Goldstone (NG) 



fields associated with ( |2.5| ) are gapless in the chiral limit and govern the low-energy physics near 



the Fermi surface. Let us label the NG modes as 

L L eSU(N f ) L /Sp(N f ) L , L R eSU(N f ) R /Sp(N f ) R , VeU(l) B , AeU(l) A . (2.6) 

For Nf = 2, Ll.r do not exist since SU(2) ~ Sp(2). First consider Nf > 4. In the presence of 
explicit breaking of chiral symmetry by the above NG modes acquire a small but nonzero 
mass. From symmetry arguments plus weak-coupling calculations, the effective chiral Lagrangian 
valid at energy scales below A is determined to be 

^=f {l^| 2 -v^|^| 2 } + ^{|5oA| 2 -v^|^A| 2 } 

f 2 r i 3A 2 r i 

+ J -j- Tr [ | d L L | 2 - vl\diL L \ 2 + (L^R)j--^ [a 2 Tr(ML R M T L[) + c.c. J . (2.7) 

The /'s are decay constants for each NG mode, and the v's are the corresponding velocities origi- 
nating from the absence of Lorentz invariance in the medium. The absence of an 0{M) term in the 
chiral Lagrangian is a consequence of the Z(2)/, x Z(2) R symmetry of the diquark pairing. 
For Nf = 2, a. similar analysis leads to 

^ = f ~2 {l^| 2 -v 2 / |^| 2 }+/ 2 ,{| )oA| 2 -v 2 ,|^A| 2 }-^{(detM)A 2 + c.c.}. (2.8) 

As V and the gluon fields decouple from the other NG modes, they will be neglected in the follow- 
ing. We note that the chiral Lagrangian at large jj, presented above is a new result. The Lagrangian 



in 1 10] is valid only for densities corresponding to jj. < m p , where m p is the p -meson mass. 



3. Partition function in a finite volume 



Next we show that the £ -regime introduced at ju = in [14] can be defined at large \x as 



well. Let us take the imaginary-time formalism and consider two-color QCD in a finite box of size 
L 4 (= V4). The masses of the NG modes due to the quark mass matrix M are denoted by wing- The 
point is that the dynamics of the system simplifies drastically in the regime 



1 r 1 

A m N G 



(3.1) 



The first inequality guarantees that contributions of non-NG modes (e.g., quarks) to the partition 



function Z are sufficiently small: e <C 1- The second inequality implies that the Compton 
wavelength of the NG modes is much larger than the linear extent of the box, which allows us 
to truncate the Hilbert space of the NG modes to its zero-momentum sector alone. Therefore the 



partition function in the £-regime (3.1) is given simply by 



Z(M)= j dA j dL L J dL R exp V 4 ^^A 2 Tr(ML R M T Ll) +c.c.} 



U(1)a SU(Nf) L /Sp{Nf) L SV(Nf) R /Sp{N f ) R 



(3.2) 
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normalized so that Z(0) = 1. The infinite-dimensional path integral is now reduced to ordinary 
integrals. For degenerate masses (M = ml), the integration can be carried out explicitly, yielding 

1 



Z(m) 



(N f - 1)!! 

where A is an Nf x Nf antisymmetric matrix with entries 



Pf(A), 



p)I p+q ( ^jV 4 A 2 m 2 



71- 



p,q 



Nf-l Nf-3 Nf-l 



(3.3) 



(3.4) 



and I p+q denotes a modified Bessel function. It is intriguing that Z in ( |3.3| ) bears strong resemblance 
to the known expression for the finite- volume partition function at /J. = in the topologically trivial 
sector [JB], (5.13)]. The latter is obtained by changing Nf and (3/7T 2 )V4A 2 m 2 in ( fO| ) and ( $A\ ) to 
2Nf and V$Lm, respectively, with £ the magnitude of the chiral condensate. It would be interesting 
to generalize our results to intermediate densities, interpolating between both results. 
For Nf = 2, explicit integration is possible for arbitrary M, with the result 



Z(M) 



-4v 4 A 2 detM ) • 

TC 2 J 



(3.5) 



4. Spectral sum rules for the Dirac operator 

In this section we review the derivation of spectral sum rules and briefly discuss their physical 
implications. Let us denote the complex eigenvalues of the Dirac operator by iX n , where the 
X n are real for /I = 0. Starting from the microscopic Lagrangian of QCD instead of the effective 
chiral Lagrangian, one may write the normalized partition function as 

Z(M) = j [9A] n'det(A 2 +M t M) e" s * / j [9A] ) % e~ s « (4.1) 



(4.2) 



measure in 



where S g = J d 4 xF^ v F^ v /4 and (• • • ) represents expectation values with respect to the 

the chiral limit. Yl' n (and later Y! n ) denotes the product (sum) over all eigenvalues with ReA„ > 0. 
We neglect the anomaly which is suppressed at large jit, and thus we assume that no zero modes 
appear. Equating ( |4.2| ) with (3^) for Nf > 4 and matching the coefficients at 0(M 2 ) and 0(M 4 ), 
we arrive at novel spectral sum rules, 

X J2 / = \ T2T2 / = \ Z TZ 



r- 



m<n '^m'^n 



0, 



r- 



X,U 4K 4 (N f -l) 2 



(V 4 A 



2\2 



(4.3) 



The vanishing of many of the spectral sums is a salient feature of the high-density limit by which 
( |Q| ) is distinguished most clearly from the conventional spectral sum rules at /I = [Q, |TJ]. We 



add that matching between (4.2) and ( pT5| ) for Nf = 2 reveals that (4.3) is correct for Nf = 2 as well. 
Introducing the spectral density p(A) and the microscopic spectral density p. s (A) defined by 



P (aW£S 2 (A-A ;i ) 



and 



p s (z) = T lim 



( TCZ. 



V 4 - 



(4.4) 
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the second sum rule in ( fO| ) is rewritten as 

[d 2 z^ = - (45) 

J z 4 4(N f -l) 2 { } 

Rez>0 

In analogy with the p = case, this formula strongly suggests that p s (z) is a universal function 
determined by the global symmetries of the problem. It leads to the observation that the smallest 
eigenvalues of @ sit at the scale of 0(1 /y/W) and that their disuibution is governed by A, in 
contrast to the situation at p = where the magnitude of the smallest eigenvalues is 0(1 /V4) and 
the quantity relevant for level correlations is not A, but £. 



Finally we comment on the generalization of the sum rules ( |4.3[ ) to the massive case. It is 
achieved by rescaling both the eigenvalues and the masses simultaneously as the volume is taken to 
infinity. We begin with the simplest case, i.e., Nf = 2 with masses m\ and wi2. In terms of rescaled 
dimensionless variables, 



Z n = K^^- , mi = m~^- , (4.6) 
K K 



the simplest sum rule reads 



1 u m 2 2 I (x) - I 2 {x) _ _ 

' • " with x = mim 2 , (4.7) 



t zl + m\// 4 /„(*) 

where ((•••)) represents expectation values with respect to the massive measure. For larger Nf the 
explicit expressions become increasingly involved. For Nf = 4 with equal masses, we have 

1 \\_2ftCy)/i(y)-3/iCy)/2Cy)+/2(y)/3(y) with y = fi? (48) 



n zl + m*// A(3h(yY-Ah{yy- + 2,h{yf) 
Similarly we can generalize p s in ( [0] ) to the double-microscopic spectral density defined by 



in terms of which we can rewrite, e.g., ( p?7| ) in the form 



nij=mj-^ — = fixed 



(4.9) 



I" d 2 z ps 2) (z;m u m 2 ) _ m 2 I Q (x) - I 2 (x) 
J z 2 + m 2 4 7 (x) 

Rez>0 1 

As emphasized in the introduction, it has been firmly established at pL = by various arguments 
[^] that the functions p s and are universal in the sense that they depend solely on the pattern 
of spontaneous symmetry breaking and not on the detailed form of the microscopic interactions. 
Hence we expect that the same holds at large p as well, even though the explicit forms of p. v and 
Ps Nf ^ defined in ( |4.4| ) and ( |4.9[ ) are still unknown. A promising approach will be to construct an 
appropriate ChRMT that corresponds to dense two-color QCD (work in progress). In this regard 
one might be tempted to consider a minimal modification of the conventional ChRMTs by addition 
of a /op-term, but this does not suffice, since in the large-p limit the latter will dominate the Dirac 
operator completely and render its dynamics trivial. 
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5. Summary 

We have constructed the low-energy effective Lagrangian for dense two-color QCD with an 
even number of flavors at large quark chemical potential ji based on the symmetry breaking pattern 
induced by the formation of a diquark condensate [[Tl]]. Also, we have identified a new finite- 
volume £ -regime for the superfluid phase at large }X. In this regime, we can exactly determine the 
quark mass dependence of the partition function from the effective theory. Matching this result 
against the two-color QCD partition function, we have derived novel spectral sum rules for inverse 
powers of the complex eigenvalues of the Dirac operator. Our sum rules explicitly show that the 
Dirac spectrum at large jj. is governed by the fermion gap A, unlike the spectrum at low ju, which 



is governed by the chiral condensate as shown in [15] 



Since the fermion sign problem is absent in this theory, our sum rules can in principle be 
checked in lattice QCD simulations. This is in contrast to real (three-color) QCD where the severe 
sign problem prevents us from observing the presumed color superconductivity, although similar 
spectral sum rules could be derived in the corresponding £ -regime Testing our sum rules for 
two colors on the lattice would enable us to measure the value of the BCS gap A at large /I for the 
first time, since previous studies of two-color QCD at nonzero /J. have only been able to determine 
the magnitude of the diquark condensate, not the gap itself. 

It would be interesting to obtain an analogue of the Banks-Casher relation at large /J. and the 
concrete form of the microscopic spectral density in the £ -regime we identified. In particular, it is 
a challenging problem to construct the corresponding random matrix model, which has turned out 
to be very successful at /J. = and small jU, to reproduce our spectral sum rules and to elucidate 
universal properties of the Dirac spectrum at large pL for both two and three colors. A detailed 
analysis of these issues is deferred to future work. 
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